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ABSTRACT 

We study winds in high redshift galaxies driven by a relativistic cosmic ray (proton) 
component in addition to the hot thermal gas component. Cosmic rays (CRs) are likely to be 
efficiently generated in supernova (SNe) shocks inside galaxies. We obtain solutions of such 
CR driven free winds in a gravitational potential of the Navarro-Frenk- White (NFW) form, 
relevant to galaxies. Cosmic rays naturally provide the extra energy and/or momentum input 
to the system, needed for a transonic wind solution in a gas with adiabatic index 7 = 5/3. We 
show that cosmic rays can effectively drive winds even when the thermal energy of the gas is 
lost due to radiative cooling. These wind solutions predict an asymptotic wind speed closely 
related to the circular velocity of the galaxy. Furthermore, the mass outflow rate per unit star 
formation rate (77^,) is predicted to be ~ 0.2 — 0.5 for massive galaxies, with masses M ~ 
10^^ — lO^^M© . We show r]yj to be inversely proportional to the square of the circular velocity. 
Magnetic fields at the /zG levels are also required in these galaxies to have a significant mass 
loss. A large t]^ for small mass galaxies implies that cosmic ray driven outflows could provide 
a strong negative feedback to the star formation in dwarf galaxies. Fiu'ther, our results will 
also have important implications to the metal enrichment of the intergalactic medium. These 
conclusions are applicable to the class of free wind models where the source region is confined 
to be within the sonic point. 
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1 INTRODUCTION 

Cosmic rays are expected to be generated and accelerated in shocks 
created by the exploding supemovae (SNe) in star forming re- 
gions. The properties of cosmic rays (CRs) in our Galaxy are well 
documented (Berezinskii et al., 1990; Schlickeiser, 2000; Dogiel, 
Schonfelder & Strong, 2002; Hillas, 2005; see also Kulsrud, 2005 
for an excellent introduction to the astrophysics of CRs). The en- 
ergy density in the proton component is about 1 eV cm~^, and CRs 
are thought to be confined to the Galactic disk for about lO'^ yr 
before escaping, presumably into the intergalactic medium (IGM). 
It is widely accepted that > 10% of the average SNe energy must 
go into accelerating the CRs so that the flux density of CRs can be 
maintained at the observed value (Hillas, 2005). 

Low energy CR protons play an important role in the ioniza- 
tion and thermal state of photo-dissociation regions in the interstel- 
lar medium of normal galaxies. As the direct interaction cross sec- 
tions are less for high energy cosmic rays, one would naively imag- 
ine that they would hardly interact with matter. However CRs in our 
neighbourhood are observed to have a very small anisotropy (about 
1 part in lO'*). This is explained by pitch-angle scattering of CRs by 
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Alfvcn waves in our interstellar medium. These waves can in turn 
be generated by the CRs themselves, by a streaming instability, as 
they gyrate around and stream along the galactic magnetic field at a 
velocity faster than the Alfven velocity (Wentzel, 1968; Kulsrud & 
Pearce, 1969; Kulsrud & Cesarsky, 1971). The waves travelling in 
the direction of the streaming CRs (and satisfying a resonance con- 
dition) are ampUfied by this streaming instability. The amplified 
waves have a wavelength of order the cosmic ray gyro-radius, and 
then resonantly scatter the CRs. Such scattering leads to transfer 
of momentum and energy from CRs to the waves. This momen- 
timi and energy is in turn transmitted to the thermal gas, when the 
waves damp out. If the CR pressure is Pc, the force exerted per unit 
volume by the CRs on the thermal gas is given by — VPc, while 
the rate of energy transfer to the gas via the intermediary of the 
waves, is given by \va ■ VPc| (Wentzel, 1971). Here va is the the 
Alfven velocity (to be defined below). The modulus is taken since 
CRs always lose energy in generating the waves, while the gas al- 
ways gains this energy as the waves damp out. In addition, power 
is transferred to the gas due to the work done by the CR pressure 
gradient in the presence of a bulk velocity v of the fluid, and this is 
given by {—v ■ VPc). 

In our Galaxy the CR pressure support plays an important role 
in the vertical force balance of the disk gas. It is also believed that 
they may regulate star formation during the evolution of galaxies 
(Ensslin et al. 2007; Socrates et al., 2008; Jubelgas et al., 2008). The 
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CR pressure gradient can also drive a wind of tliermal gas and CRs 
from tiie galaxy. The possibility that CRs play an important role in 
driving outflows from our Galaxy has been widely studied in the 
literature (Ipavich, 1975; Breitschwerdt et al., 1987; 1991; 1993; 
Zirakashvili et al., 1996; Ptuskin et al., 1997; Breitschwerdt et al., 
2002; Everett et al., 2008; 2009). Such models can also explain bet- 
ter the observed diffuse soft X-ray emission of our Galaxy (Everett 
et al., 2008). High redshift galaxies are expected to have a much 
larger rate of star formation than our Galaxy. Therefore, it is natu- 
ral to consider the possibility of their developing outflows driven by 
cosmic rays. This forms the main motivation for the present work. 

Furthermore, observations of high-z Lyman break galaxies 
(LBGs) frequently show evidence for galactic scale superwinds 
(Pettini et al., 2001). Profiles of C IV and O VI absorption lines 
seen in the high redshift damped Lyman-a systems (DLAs) are 
also consistent with them originating from outflows in DLA pro- 
togalaxies (Fox et al., 2007a, 2007b). The possible importance of 
outflows in the high redshift universe is also indicated by the pres- 
ence of metals in the IGM as traced by the Lyman-a forest (Tytler 
et al., 1995; Songaila & Cowie, 1996; CarsweU et al., 2002; Sim- 
coe et al., 2002; Bergeron et al., 2002; Ranch et al., 1997; Songaila, 
2001; Ryan- Weber et al., 2006). These metals can only originate in 
stars and have to be then transported out into the IGM via outflows. 
Thus outflows appear to be ubiquitous in high redshift galaxies. 

There exist several semi-analytical and numerical models of 
outflows aimed at explaining the presence of metals in the IGM 
(Nath & Trentham, 1997; Efstathiou, 2000; Madau, Ferrara & 
Rees, 2001; Furlanetto & Loeb, 2003; Scannapieco, 2005; Samui, 
Subramanian & Srianand, 2008). All these models either use the 
pressure or momentum of the thermal gas, generated from SNe ex- 
plosions to drive the outflow. There have also been models based 
on radiation pressure driving the winds (Dijkstra & Loeb, 2008; 
Nath & Silk, 2009). The role played by cosmic rays (which will be 
copiously produced in such high redshift starburst galaxies) has not 
been explored so far in the context of high redshift galactic winds. 

Moreover, it is well known that for a 7 = 5/3 thermal gas, 
transonic wind solutions with the fluid velocity being accelerated 
from subsonic to supersonic speeds exist only if there is some form 
of energy and/or momentum input till the critical point (Lamers 
& Cassinelli, 1999). An example of such a solution is that due to 
Chevalier and Clegg (1985), where one assumes that mass and ther- 
mal energy are injected into the flow upto the critical point. The 
Chevalier and Clegg (1985) solution also ignores gravity, assuming 
the gas is hot enough not to feel the gravitational influence of the 
galaxy. It appears unnecessarily restrictive to assume that such a 
situation of energy input till the sonic point can always be obtained 
in a galaxy. On the other hand cosmic ray driving can naturally 
provide the required energy and/or momentum input, even if the 
energy sources (i.e. SNe) are well within the critical point, or the 
thermal energy is lost due to radiative cooling. This forms an ad- 
ditional motivation for considering such cosmic ray driven winds 
from high redshift star forming galaxies. 

In the present work, we construct models of CR driven winds 
from star forming galaxies. Although our emphasis is on high red- 
shift galaxies, our wind solutions would be equally applicable to 
any galaxy where the star formation results in copious production 
of CRs. The pioneering work by Ipavich (1975) considered cos- 
mic ray driven wind solutions in a point mass gravitational poten- 
tial. Here, we closely follow the methodology of Ipavich (1975). 
However we take the gravitational potential to be described by the 
standard NEW profile (Navarro, Frenk & White, 1997), relevant 
for galaxies. The paper is organized as follows. In the next sec- 



tion we recall the equations for steady, spherically symmetric, cos- 
mic ray driven outflows, and also derive the general wind equa- 
tion. This allows us to estimate the mass loss due to such out- 
flows. Section 3 casts these equations in dimensionless form suit- 
able for numerical solution. A suite of such solutions of cosmic ray 
driven outflows from high redshift galaxies with masses ranging 
from 10^ — 10^^ Mq, are presented in section 4. For the smaller 
mass galaxies, it turns out that the gas may radiatively cool signif- 
icantly, and so in section 5 we consider a suite of cold, cosmic ray 
driven, wind solutions. We end with a discussion of our results in 
section 6. 



2 THE CR DRIVEN OUTFLOW 

Let us consider a steady spherically symmetric outflow of thermal 
gas of density p and 7 = 5/3, and ultra relativistic cosmic ray 
particles, with 7 = 4/3. We assume that the cosmic ray component 
can be described as a fluid with negligible mass density, but having 
appreciable energy density Uc and a pressure Pc = IJcl'i- The 
spherically symmetric form of the fluid and cosmic ray equations 
are given by (cf. Ipavich, 1975), 



pvr = q — constant, 



dv 
pv— — 



dr 



1 d 

J.2 



pvr 



dP 
dr 

5 P 



dPc 
dr 



GM{r)p 



1 2^ 

-V H 

2 2 p 



-pvr 



,GM{r) 



+ 1, 



1 A [4Py{v + VA)] =~I. 
r^ dr ^ 'J 



(1) 



(2) 



(3) 



(4) 



Here v and P are the radial velocity and pressure of the thermal gas, 
VA ~ B I (A,-K pY^"^ is the Alfven velocity, with B describing the 
radial component of the magnetic field. As in Ipavich (1975), we 
assume the magnetic flux, given by B{r) (47rr^) to be constant. The 
equation of continuity describing mass conservation, is integrated 
to give Eq. [T| The momentum equation (Eq. |2j, includes pressure 
gradients due to both the thermal and cosmic ray components, and 
the gravitational attraction due to a mass distribution M (r). Energy 
conservation is described by Eq.[3] where the / term. 



/ = - (t; + va) 



dr 
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incorporates the exchange of energy between the cosmic ray and 
thermal gas components. The cosmic ray evolution itself is de- 
scribed by Eq.|4] 

We can also add Eqs. [3] and |4] describing respectively energy 
conservation and CR evolution, and integrate the resulting equation 
to get another algebraic conservation law. 
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where Ce is some constant of integration. If we multiply the above 
equation by 47r, then the terms on the LHS are respectively the 
rate at which, the kinetic energy, enthalpy of the thermal gas, grav- 
itational energy and cosmic ray 'enthalpy' are flowing through a 
sphere of radius r. The sum is a constant at any radius and equal 
to 4nCe. In the inner regions of the wind, the kinetic energy con- 
tribution will be small, and the thermal and CR pressures contri- 
butions will have to be sufficiently large to overcome the gravity. 
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As r oo, all terms except the kinetic energy contribution will 
go to zero, and Ce q{y'i^/2), where Vac is the wind velocity as 
r — > oo. 

A straightforward manipulation of the above equations, which 
we describe in Appendix|A| allows us to derive the wind equation. 



2v 



1 



GM(r) 



2rci 



Av 
dr 



where the effective sound speed c, is determined by. 
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One sees that flow has to pass through a critical point where both 
the numerator and denominator of Eq.|7]has to vanish. Indeed the 
numerator has to go to zero faster than the denominator, to obtain a 
regular solution for the flow. At the critical point, r — I, the velocity 
of the gas is given by. 



v{l) = c.(0 = y/GM{l)/2l = Vc 



(9) 



where we have defined the circular velocity of the galactic potential 
at the critical point r = ItobeVc ^ (GM (I) jlY''^ . Note that Vc 
is also roughly the circular velocity at the virial radius r = r„i,,, 
if the galaxy has approximately a flat rotation curve; that is 14 ~ 

We can now estimate the mass outflow rate as follows. Let the 
total luminosity that is being pumped by the SNe to the outflowing 
material be Lq, and assume that only a fraction / of this ends up 
as the wind kinetic energy. This fraction depends on the gravita- 
tional potential of the galaxy and the distribution of mass and en- 
ergy sources with respect to that potential. In Appendix|B]we have 
given an estimate of / for a uniform source distribution and find a 
typical value of / ~ 0.1. Moreover by solving the fluid equations 
numerically, we will generally find that the asymptotic velocity of 
the outflowing gas as r — > oo, «oo, is closely related to the velocity 
at critical point and hence the circular velocity, i.e. 
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with F/ \/2 is of order unity. Then, 

47ra = ^M„t)L = fLo = f{vMsF){e^EsNe). (11) 

Here ~ 4:nq is the mass outflow rate, Msf is the star forma- 
tion rate, v is the number of SNe produced per unit mass of stars 
formed, EsNe = 10^^ erg is the energy produced per SNe and 
is the fraction of that energy that drives the outflow]^ Hence the 
mass outflow rate is given by. 
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The factor depends on the initial mass function of the stars 
formed, and is ~ 1/130 Mq ^ for a Salpeter IMF with a lower and 
upper cut-off masses of miom =0.1 M© and m^p = 100 Mq. 



Changing miow to 1 Mq gives v = 1/50 Mq^. We will adopt 
e„ ~ 0.1. The wind mass loading factor is then, 
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where for the numerical estimate, we have scaled by some typi- 
cal values of e^, v and Vc- Our estimate in Eq. ^] leads to one 
of the important results of our work, that the mass loading factor 
is expected to increase with a decreasing galaxy circular velocity. 
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V- 



w^,,. For a massive galaxy with 
Vc like our Galaxy, and taking / ~ 0.1, F ~ \/2 (see below), 
we expect to get ri^ ~ 0.2. On the other hand CR driven outflows 
from dwarf galaxies with much smaller Vc are expected to lead to 
a much higher mass loading factor. The numerical solutions which 
are obtained below, bear out these general expectations. 



3 NUMERICAL SOLUTIONS 

We now consider numerical solutions of the fluid equations includ- 
ing the CR component. Note that during most of its evolution the 
outflow traverses through the dark matter halo of a galaxy, whose 
potential is generally well described by the standard NFW form 
(Navarro, Frenk & White, 1997). The mass distribution associated 
with the NFW potential is given by, 



where the function 5(2) is given by, 
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Here M is the total galaxy mass, r^ir the virial radius of the halo, 
and c the concentration parameter, which decides the radius be- 
yond which the rotation velocity transits to an almost flat form. 
We choose c to be larger than what would obtain for just the dark 
halo, in order to mimic the fact the observed rotation curves of 
disk galaxies are flat well within their optical radius. Therefore, our 
model NFW potentials provide a reasonable approximation to the 
gravitational potential of the galaxy, at least in the spherical limit. 

In order to obtain the wind solutions,we first transform all the 
equations to a dimensionless form. We use the location of the crit- 
ical point, say I, as the length variable and define a dimensionless 
radius x = r/l. Thus the critical point occurs at 2: = 1. Using 
the constants q, Cc and I we define the dimensionless density (p*), 
velocity (u), thermal pressure (Op) and cosmic ray pressure (Sn) 
respectively as 
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In terms of these dimensionless variables, the evolution equa- 
tions |TJ, |2|, |6| and |3| respectively, for the CR driven outflow 
become 



* 2 o 
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^ A qualitatively similar estimate, without the factor /, is given by Murray, 
Quataert & Thompson (2005). 
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Here we have also defined the dimensionless parameters of the sys- 
tem, 
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The parameter 'a' represents the ratio of a typical gravitational po- 
tential energy to thermal plus CR energy, while the parameter b 
gives a measure of the strength of the magnetic field. These defini- 
tions are similar to that adopted by Ipavich (1975). The parameter 
c' gives the dimensionless measure of the location of the critical 
point, in terms of the scale length r^i^jc associated with a dark 
matter halo. This parameter was not present in the work of Ipavich 
(1975), as there gravity is due to a point mass. Its presence increases 
the complexity of finding transonic solutions. Note that we need to 
solve only two differential equations, due to the availability of the 
two algebraic constraints representing mass and energy conserva- 
tions. We first solve for gas pressure, using Eq.|19[ 
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Substituting this in Eq.|20| we obtain 
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Next we substitute Eq.|24]and Eq.|25]into Eq.[T8]to get the velocity 
gradient: 
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There is again a critical point to the solution where both the numer- 
ator and denomination in Eq. |26] vanish. These conditions imply 
that at the critical point, where a; = 1, the velocity u = Uc, where 
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The cosmic ray pressure at the critical point, 6fi, can be obtained 
using Eq.|27|and|29[ 

" \ 1 + c' J [6u'i+25u,{ucb)^/^ + Uucb]' ^ ' 

Since both numerator and denominator vanish at the critical point, 
we use L'Hospital rule to evaluate the velocity gradient (i.e. Eq.|26|l 
there. This results in a quadratic equation for du/dx at the critical 
point given by, 
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with all quantities in the above three equations evaluated at a; = 1 
and u = Uc - Note that A = 2ucB. 

Our numerical solutions start from the critical point which is 
at a; — 1. At the critical point both the numerator and denomina- 



tor of Eq. 26 vanish though du/dx is finite. Hence, we calculate 
du/dx at the critical point from Eq.jsT] after adopting specific val- 
ues for the parameters a, b and c'. We take the positive root of this 
equation since we want a subsonic flow to transit to a supersonic 
flow. We go both forward and backward in x and solve the coupled 
differential equations, Eq.|25]and Eq.|26] This gives the dimension- 
less velocity u{x) and the dimensionless CR pressure 6n{x). Then 
we use the algebraic Eqs. [TT] and |24] to obtain the dimensionless 
p, (x) and 9p{x), respectively. These dimensionless quantities can 
be converted to the dimensional quantities (p, v, P, Pc), using the 
transformations of Eq.[T6| 



4 RESULTS 

The numerically computed CR wind solutions for several different 
galactic masses, 10* Mq to 10^^ Mq, and different parameters, are 
given below. Table [T] summarizes the assumed parameters and the 
derived physical quantities associated with each of these solutions. 
The double solid line in the table separates the assumed parameters 
(above this line) and the derived quantities (below the line). 



4.1 CR driven wind from a massive galaxy 

First consider the example of wind evolution in the potential of a 
10^^ Mq galaxy described by the NEW density profile. We model 
the parameters of this potential to mimic outflows from a massive 
high redshift galaxy, by adopting c = 15 and a virial radius r^ir ~ 
100 kpc. We take the critical point to occur at I = r^ir/W or 
I — 10 kpc. This fixes the value of c = 1.5. We also take a = 2.0 
and b = 1.5 and solve the dimensionless wind evolution equations. 
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In order to convert the wind solution to a dimensional form, 
we also have to adopt a value for Ce. Recall that inCe = fLo, 
with / ~ 0.1 (Eq.[TT|and Appendix|B](. Here, as mentioned earlier, 
Lo is the rate of energy input by SNe in the galaxy, and is given by, 

^o = 3xW [yj,j^)[^i)[-^)^^SS .(34) 

For the solution shown in Fig.[T| we have assumed a star formation 
rate (SFR) Msf ~ (10//) MQjr''^, which leads to 4nCe = 
3 X 10*^ erg s~^. This SFR is typical of that obtained in massive 
high redshift Lyman break galaxies (Erb et al., 2006). The assumed 
and derived parameters are summarized as model A in Table[T] 

The fluid velocity obtained for this model is shown in panel (a) 
of Fig.[T]as a solid line. The wind starts from a finite subsonic veloc- 
ity (v < c) at ro — 7.8 kpc. It passes through the critical point at 
r = 10 kpc and the velocity continues to increase before reaching 
an asymptotic velocity ~ 292 km s^^. This should be compared 
with the circular velocity at the critical point, Vc ~ 270 km s^^ 
and thus for this solution, F/%/2 ~ 1.07. Panel (a) of Fig.[I]also 
shows the circular velocity of the galaxy (as a dotted line), and 
the Alfven velocity (va) (as a dashed line). The circular velocity 
is reasonably constant, varying by about 10% between 5 and 50 
kpc. From the figure it is also clear that the Alfven velocity de- 
creases with radius after a small initial increase (like in the solu- 
tions given by Ipavich (1975)). The cosmic ray pressure and the 
thermal pressure are shown by solid and dashed line respectively 



in panel (b) of Fig. [ij At the critical point Pc ~ 0.50 eV cm~ 
and P ~ 0.59 eV cm~'^. In the present solution, the CR pressure is 
comparable to the thermal gas pressure at the base of the flow, then 
goes below the thermal pressure and then again becomes higher at 
large radii. This latter feature arises because 7 = 4/3 for the cos- 
mic ray component whereas 7 — 5/3 for the thermal gas. As the 
wind expands out, the thermal gas pressure then decreases faster 
with radius compared to the cosmic ray pressure. 

In panel (c) of Fig.[T]we show the density of the wind material. 
At the base of the outflow the density is ~ 0.008 cm~^ and falls 
of asymptotically as r^^, once the velocity of the flow becomes 
nearly constant. Note that if we have a mass of lO^^A/0 of baryons 
in a galaxy distributed within 10 kpc, the average density would be 
1 cm""'. Thus the gas density required in the wind solution is mod- 
est in comparison. The resulting temperature at the critical point is 
2.6 X lO'' K and falls of monotonically with radius thereafter. 

Finally, we show the magnetic field strength in the wind ma- 
terial in panel (d) of Fig.[T| Recall that we assumed B{r) oc r^^. 
The proportionality constant is fixed from the parameter b of the 
model. We find the value of B at the critical point to be ~ 6.2/iG. 
This is comparable to typical values of a few to lO/iG, obtained 
for the mean field in nearby spiral galaxies (Beck, 2008), whereas 
it is smaller than the total fields of 50 — 100/iG inferred in some 
nearby starburst galaxies (Chyzy & Beck, 2004). We find that a 
threshold minimum value of the magnetic field is required to ob- 
tain a cosmic ray driven wind passing through the critical point. If 
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the field were smaller than this threshold, the cosmic ray pressure 
becomes negative at some point. For the dimensionless parameters 
a = 2.0, c' = 1.5 that we have adopted, we require the parameter 
determining the magnetic field strength b > 0.5 to obtain a tran- 
sonic wind solution. Thus one still requires /iG fields to have CR 
driven transonic winds. This requires efficient dynamo generation 
mechanisms to operate in high redshift galaxies (cf. Brandenburg 
& Subramanian, 2005; Shukurov, 2007; Sur, Shukurov & Subra- 
manian, 2007). Recent observations of Faraday rotation towards 
quasars showing Mg II absorption at high redshift z ~ 1 — 2, 
suggests that organized magnetic fields of high strengths ~ 10^ 
G, could be associated with normal galaxies even at these redshifts 
(Bemet et al., 2008; Kronberg et al., 2008). 

We should also mention in passing that the magnetic field need 
not be coherent on the galactic scales, for the CRs to transfer their 
energy and momentum to the thermal gas. However one does need 
some degree of coherence. First, for Alfven waves to be resonantly 
amplified, the field should be coherent over scales much larger than 
the CR gyro radii (r-g), which is of order lO^'^e(GeV) cm, where e 
is the CR energy in units of GeV (Kulsrud, 2005). This scale is 
clearly much smaller than typical scales in the galactic wind. Sec- 
ond, energy transfer via the waves occurs only if either the forward 
or the backward travelling wave (with respect to the field direction) 
dominates the wave spectrum (Skilling, 1975, Eq. 9). Note that the 
Alfven waves basically travel along the field lines as they damp 
out, and unless this damping length is smaller than the coherence 
scale of the magnetic field, waves going in one direction would 
meet oppositely travelling waves, and the net wave heating would 
be negligible. For the physical conditions in the outflow solutions, 
we find this damping length is generally smaller than a fraction of 
a parsec (cf. Kulsrud, 2005, Chapter 12). Thus we need magnetic 
field coherence scales to be larger than this value. We also note that 
momentum transfer (without the corresponding energy transfer via 
waves) imposes a milder constraint that the field coherence scale be 
only much larger than rg. Note that radio observations of several 
edge on galaxies do show significantly coherent radial fields ema- 
nating from the galaxy (Chyzy et al., 2006; Heesen et al, 2009a, b; 
Krause, 2008). 

The main advantage of our work, is that one can predict a mass 
outflow rate given other parameters. The restriction on the mass 
outflow rate arises because the solution is constrained to be regular 
at the wind critical point. In this particular example of a 10^^ Mq 
galaxy, we find the mass outflow rate is Anq — 16.2 Mq yr~^. 
This should be compared with the required star formation rate of 
MsF ~ (10//)MQyr-\ giving n^u ~ 1.6/ ~ 0.16 [for / ~ 
0.1]. 

Note that the mass outflow rate is directly proportional to the 
value of Ce that one assumes. A value of Ce larger by a factor 
10, would also lead to 10 times higher value of g, for the same 
value of a. This can be directly seen from Eq. |^ The value of 
rju, however remains the same as both q and AisF increase by the 
same factor. Our model B is an example of such solution where we 
assume 10 times higher inCe = 3 x 10^^ erg s~^, keeping all 
other parameters same as in model A. Increasing Ce by a factor 10 
leads to 10 times higher values of P, Pc and p and B higher by a 
factor 



4.2 Wind solutions for different galaxy masses 

We now consider wind solutions for a range of galaxy masses 
adopting the same dimensionless parameters, a = 2.0, b = 1.5 
and c' = 1.5. For each galaxy mass we fix a virial radius assuming 



the galaxy (and its dark matter halo), collapsed from a 3cr density 
fluctuation in a standard LCDM cosmology. The critical point is 
taken to be at Z = rui^/lO. Further we assume that the energy out- 
flow rate Ce is proportional to the galaxy mass, and is given by 
AnCe = 3 X 1O''^(M/1O^^M0) erg s"\ Such a scaling assumes 
that a fixed fraction of the galactic mass is converted to stars over 
a fixed timescale. Thus the star formation rate and hence the rate 
of energy input by SNe is proportional to the mass of the galaxy 
(normalised as earlier to a SFR of ^ (10//)Mq yr~^ for a mas- 
sive lO^'^M© galaxy). The derived physical parameters of the wind 
solutions, for galaxy masses ranging from 10^^ — 10* Mq are pre- 
sented in Table[T] as models C to F respectively. 

Several trends in the derived physical parameters can be seen 
from Table[T] As the galaxy mass decreases from 10^^ — 10* Mq, 
the gas and cosmic ray pressures at the critical point increase by a 
factor ~ 4, from about 0.5 — 2 eV cm~^, while B increases by a 
factor ~ 2 from 6.2 — 13.2 /iG. There is a much larger change in the 
density of the wind material, which increases by nearly a factor ^ 
10^ from ~ 0.003 - 2.1 cm as the galaxy mass decreases by 4 
orders of magnitude. Correspondingly the temperature of the wind 
which is proportional to P/ p, decreases from about 2.6 x 10^ K, 
for a 10^^ A/0 galaxy to T ~ 1.5 x 10* K for a 10* Mq galaxy. 
The increased density and decreasing temperature associated with 
these winds for lower mass galaxies implies that gas cooling would 
become more and more important for them. We will return to this 
aspect below. 

Further, for the models A to F in Table[T| we obtain very sim- 
ilar values of F, with F/y/2 ~ 1.05 - 1.08. Therefore, in all these 
cases the asymptotic wind speed follows closely the circular veloc- 
ity of the galaxy at the critical point, with the factor F/ \/2 of order 
unity. There can however be a factor of 2 scatter in this relation for 
the range of dimensionless parameters a, b and c' that we consider 
(see models E, G, H and I of Table[TJ. Moreover, the mass loading 
factor follows very closely the scaling with circular velocity at the 
critical point predicted by analytic arguments, 77m// oc V^^ . Thus 
while massive high redshift galaxies have a mass loading factor of 
^ 0.2, small mass galaxies, with total mass M ~ 10'' - 10* Mq 
have a large rj^ ~ (7.6 — 28) for / ~ 0.1. In Samui, Subramanian 
& Srianand (2008, 2009) we found that the volume filling of the 
IGM by outflows is mainly dominated by small mass galaxies with 
M < 10^ Mq. Hence the results on small mass galaxies obtained 
here have very important implications for such models. 



4.3 Varying the dimensionless parameters 

Let us now consider the effect of dimensionless parameters a, b and 
c'. To begin with, we concentrate on their impact on the dimension- 
less velocity profile. In all transonic solutions u starts from some 
finite value uq and at a finite x — xq and increases till it passes 
through the critical point at a; = 1. Note that the mass conservation 
equation, Eq. [TT] implies that the solutions have to start at finite 
X = xo, and with a finite velocity, to avoid infinite densities. This 
feature can be changed if we have source terms in both the mass 
and energy conservation equations at least upto x = xq. One of 
our aims is also to look for transonic solutions where xo is as small 
as possible, such that sources need not extend to the critical point 
x = l. 

we keep a fixed b = 1.5 and c' — 1.5 and vary a. 



In Fig. 



4.1 



The starting point xo decreases as one adopts larger and larger a, 
till about a — 2, where it reaches xo — 0.78. Further increase in a 
does not lead to a decrease in xq. Also as can be seen from Fig. |4.1| 
Uc increases with a, a feature which follows from Eq. [28] which 
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Figure 2. The dependence of dimensionless velocity u on the dimensionless model parameters a, b and c' are shown in panels (a), (b) and (c) respectively. 
Values of parameters that remained fixed are given on the top of each panel. Vertical dotted lines in each panel represent the critical point x = 1. 



Table 1. The summary of various wind models. 











Models 










Parameters 


A 


B 


C 


D 


E 


F 


G 


H 


I 


a 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


0.5 


2.0 


2.0 


b 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


4.0 


1.5 


c' 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


2.5 


mass (Mq) 


1012 


1012 


10" 


10"' 


109 


10« 


109 


lO'J 


109 


r^ir (kpc) 


100 


100 


30.7 


11.9 


4.7 


1.7 


4.7 


4.7 


4.7 


/ (kpc)t 


10 


10 


3.1 


1.19 


0.47 


0.17 


0.47 


0.47 


0.78 


c 


15 


15 


15 


15 


15 


15 


15 


15 


15 


iirCe (erg s^^) 


3 X 10*1 


3 X 10^2 


3 X lO-io 


3 X 10^9 


3 X 10=** 


3 X 10^7 


3 X 10=*** 


3 X 10^** 


3 X lO^s 


ro (kpc) 


7.8 


7.8 


2.4 


0.92 


0.36 


0.13 


0.43 


0.30 


0.35 


P(eVcm-3) 


0.59 


5.9 


1.10 


1.44 


1.84 


2.66 


0.55 


2.19 


1.09 


Pc (eV cm" 3) 


0.50 


5.0 


0.93 


1.22 


1.56 


2.25 


0.46 


1.08 


0.71 


BipG) 


6.2 


19.8 


8.5 


9.8 


11.0 


13.2 


7.8 


18.0 


7.5 


n (cm"'^) 


2.7 X 10"^ 


0.027 


0.015 


0.077 


0.388 


2.03 


0.097 


0.388 


0.231 


T(K) 


2.6 X lO'' 


2.6 X 10^ 


8.5 X 10^ 


2.2 X 10^ 


5.5 X lO-i 


1.5 X 10"' 


6.5 X iC 


6.5 X 10'* 


5.4 X 10* 


Vc (kms-l)t 


271 


271 


155 


79 


40 


21 


40 


40 


40 


Voo (km s^i) 


292 


292 


166 


85 


42 


22 


66 


38 


31 


F/V2 


1.08 


1.08 


1.07 


1.08 


1.05 


1.05 


1.65 


0.95 


0.78 


iirq (Mq yr-i) 


16.2 


162 


4.99 


1.93 


0.76 


0.28 


0.19 


0.76 


1.27 




0.16 


0.16 


0.5 


1.9 


7.6 


28 


1.9 


7.6 


12.7 



t I = r^ir/lO 

t calculated at r = ( 

§ rjnj is calculated taking / = 0.1 
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predicts Uc oc a^''^. Importantly, we find that transonic solutions 
are not possible for arbitrarily small values of a. For example one 
requires a > 0.5 to have transonic solutions adopting 6 = 1.5 and 
c' = 1.5. 

The effect on the velocity profile of varying h and c are shown 
respectively in Fig. 4.1 and Fig. 4.1 As one increases b or c keep- 



quiring a larger magnetic field, but starts the flow at smaller radius 
compared to the critical point. 



ing other parameters fixed, xq decreases, and the starting point for 
the wind can be at a smaller radius. Thus a larger value for the mag- 
netic field or a more concentrated galaxy potential aids in starting 
the outflow at a smaller radius At the same time, a minimum b 
and c' are required to get a transonic solution, with other parame- 
ters fixed. 

We now consider the changes to the wind solutions in dimen- 
sional form, which obtain on changing the adopted dimensionless 
parameters. Our models G, H and I study this aspect. We focus on 
the case when the galaxy mass is 10^ Mq because of the impor- 
tance of such small mass galaxies to the volume filling of the IGM. 
We find that decreasing the parameter a by a factor 4 (i.e. model G) 
basically decreases the mass outflow rate by the same factor. This is 
as expected from Eq.|2T] where we see that g oc a. The decrease in 
a also brings the base of the outflow, ro, closer to the critical point 
I. Model H shows that increasing b basically leads to a solution re- 



^ Such changes in zq also alters /, but from our work in Appendix [b] by 
at most a factor of ~ 1.3. 



Finally, we have illustrated the effect of increasing the param- 
eter c' from 1.5 to 2.5 in model I. In Fig. [s] we have shown the 
characteristics of this model. Since we wish to keep the NFW po- 
tential form similar, we have adopted the same concentration pa- 
rameter c, which implies that the critical point has to occur at a 
larger radius; l/r^ir = 0.17 for this model instead of 0.1. As also 
mentioned earlier, we see that increasing c' allows the wind to start 
at an even smaller tq/I ~ 0.45. It also leads to a larger mass out- 
flow rate, with Aixq = 1.28M0 yr~^, larger by a factor of 1.7 
compared to the solution with c' = 1.5. The cosmic ray pressure 
again begins slightly smaller than the gas pressure at the base of the 
outflow r = ro = 0.35 kpc, with values of Pc = 0.7 eV cm~"^ and 
P = 1.1 eV cm~^ at the critical point r = 0.78 kpc, but exceeds 
the gas pressure asymptotically. The asymptotic wind velocity is 
31 km s^^, and in this case smaller than the circular velocity at 
the critical point, but comparable to the asymptotic circular veloc- 
ity. The wind density is ~ 0.23 cm~'^ at the critical point and then 
falls as 1 /r^ asymptotically. The magnetic field strength is of order 
7.5^G at the critical point. 
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Figure 4. Temperature profile for the models A and E of Table[T]are shown 
in panel (a) and (b) respectively. 



4.4 Cooling of the wind material 

The wind solutions presented above do not take into account ra- 
diative cooling of the thermal gas. Note that in these solutions 
the density of the wind material at the critical point is n{l) ~ 
0.003 — 2 cm^'^ for galaxy masses in the range 10^^ — 10* Mq. 
The corresponding temperatures of the wind at the critical point are 
2.6 X lO'' — 1.5 X 10* K for this mass range of galaxies. We show 
the temperature profile for the model A and E of Table [T] in top 
and bottom panels of Fig. [4] respectively. The radiative cooling of 
the gas depends on the temperature, density as well as the metal- 
licity of the gas. Taking a typical value for the cooling function 
Ajv = 5 X 10"^^ erg cm^ s"^ (Sutherland & Dopita, 1993), the 
cooling times are tc ~ 10* yr and tc ~ 10'' yr for 10^^ Mq and 
10* Mq galaxies respectively. This should be compared with the 
advective time scale of the wind at the critical point, ta ~ l/v{l). 
We find ta ~ 5 X 10^ yr and 10* yr for the lO" Mq and 10* Mq 
galaxies respectively. Therefore, radiative cooling is very important 
especially for low mass galaxies. However the cosmic ray pressure 
by itself could drive a wind even if the thermal gas cools rapidly, as 
was found by Ipavich (1975), for outflows in a point mass potential. 
We therefore examine in the next section, if one can obtain similar 
cold wind solutions (where the thermal pressure is neglected), from 
galaxies in a NFW halo potential. 



5 COLD WIND SOLUTIONS 
5.1 Basic formalism 

We assume that radiative cooling is so efficient that most of the 
thermal energy is lost from the gas, and so neglect the influence 
of the thermal pressure compared to the CR pressure in the mo- 
mentum equation. However, cosmic rays can still transfer momen- 
tum (via Alfven waves) to the thermal gas. The equations of mass 
conservation (i.e. Eq.[T|(, momentum conservation (Eq.[2]( with the 
thermal pressure neglected, and the equation governing cosmic ray 
evolution (Eq.[4]and Eq.|5j still hold. And they form by themselves 
a closed system describing the evolution of p, v and Pc- Again we 
assume that the magnetic field falls of as 1/r^. From Eq. [T| Eq. |4] 
and Eq. [3] we still get Eq. |A8| of Appendix A. Substituting this in 
Eq.[2] neglecting P compared to Pc, gives 



dv ,2 dp 

pv- he, — = 



dr 



dr 



pGMjr) 



where now 



4 (i; + va/2) 



3 p {v + va) 



1/2 



(35) 



(36) 



is the effective sound speed. Eliminating dp/dr using Eq. 
obtain the cold wind equation 



A7 



dv 
dr 







GM{r) 


2v 


1 


2rc'? 


r 

















(37) 



Hence, the cold wind equation (Eq. [37} has the same form as the 
hot wind equation (Eq.|7|; only the effective sound speed has now 
changed since we are neglecting the thermal pressure of the gas. 
Thus the arguments relating the mass loading factor and the circular 
velocity of the galaxy given in section 2, and leading to Eq. [13] 
are still valid, except that effective / is likely to be smaller due to 
radiative cooling. 

In order to obtain numerical solutions, we again use the trans- 
formations of Eq.[T6] and get 



p ux 



2, 



2— + — 
dx ux'^ 



ln(l + cx) 



1 + dx 



+ x —— = 0, 
dx 



d^ 
dx 



u + 



[uh) 



1/2 



U + 



(ub) 



1/2 



dBn 
dx 



(38) 



(39) 



(40) 



Basically Eqs. |38| and[39|come from Eqs.[T]and[2]and Eq.|40| comes 
from combining Eqs.[4]and[5] Note that in this case Ce used in the 
transformation equations (Eq.|16|) is just an arbitrary constant and 
does not represent the conserved luminosity of the system like in 
Eq.[6| 

Following the same procedure as in section 3, we arrive at two 
simultaneous differential equations for the flow: 



dBn 
dx 



4 


"2 Idu 

X udx 




[ 2. J 




u -\ 

X 





(41) 
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Figure 5. Characteristics of cosmic ray driven cold wind coming from a 10^^ Mq galaxy (thick lines). We show velocity, thermal gas pressure, cosmic ray 
pressure, the density of the thermal gas and the magnetic field as a function of radial distance from the galactic centre. The vertical dotted lines in each panel 
represent the critical point. For comparison we also show the corresponding hot wind solution with thin lines. 



and 



du _2u ^\ 1 + c'x 

dx X 



where 



f — -Ouux^ 



{ub) 



1/2 



2x 



(ub) 



1/2 



(42) 



(43) 



Note that Eq.|4T|and Eq. [25]are identical (a feature also found by 
Ipavich, 1975). 

The flow again has a critical point when both numerator and 
denominator of Eq.|4T]vanish and we again choose the critical point 
to be at X = 1. At the critical point, 



ln(l + c') - 



1 + c' 



and 



u + {uby^^ 



u + 



(ub) 



1/2 



(44) 



(45) 



Using the L'Hospital's rule, we evaluate the dimensionless velocity 



gradient at the critical point. 



du V 



du 



(46) 



where 



A = 10/3^31/3 + 28, Z3 = 4^(4-5/3-2/3^), 

C = -4u^(5/3^ + 14/3 + 2) + f (1 + /3)(2 + /3), 

(1 + c')^ 

13 = va/v = {ub)^^"^ /u, 

with u in Eqs.|45]and|46]is evaluated at x = 1. We now turn to the 
numerical solution of the cold wind equations for a range of galaxy 
masses. 



5.2 Cold wind solution from a massive galaxy 

Consider first CR driven cold wind from a massive, 
galaxy. We choose all the parameters, except Ce, as in model A. 
Recall that the value of Ce is an arbitrary constant for the cold wind 
solutions. We fix its value by demanding that at the critical point, 
r = I = IQ kpc, the cosmic ray pressure is similar to that which 
obtains for the hot wind solution described in the previous section 
(see model A of Table[TJ. We do this, so as to directly compare the 
mass outflow rate (47rg) in the two cases. For this purpose, it turns 
out that one needs to adopt a value for 47rCe = 1.5 x 10^^ erg 
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Figure 6. Characteristics of cosmic ray driven cold wind coming from a 10^ Mq galaxy (thick lines). We show velocity, thermal gas pressure, cosmic ray 
pressure, the density of the thermal gas and the magnetic field as a function of radial distance from the galactic centre. The vertical dotted lines in each panel 
represent the critical point. For comparison we also show the corresponding hot wind solution with thin lines. 



half the value used in the hot wind solution of Fig.fT] The resultant 
dimensional cold wind solution is shown in Fig. pTas thick lines. 
We also show for comparison, the corresponding hot wind solution 
given in Fig.[T] as thin lines. 

The first feature to note (panel (a)) is that, the cold wind solu- 
tion can start from an arbitrarily small radius and velocity; however 
from mass conservation the density would then blow up. Having an 
upper limit on the density will then decide the location of the base. 
Such a feature is also seen in the cold wind solutions obtained by 
Ipavich (1975). This is unlike the hot wind solution, where the flow 
starts from a finite radius with velocity significant fraction of the 
local circular velocity. The final asymptotic velocity of the cold 
wind solution turns out to be higher than the corresponding hot 
flow, but still close to the circular velocity at the critical point, with 
= 1.26. The cosmic ray pressure through out the flow is 
similar to that of previous solution. However, both the density and 
the magnetic field are smaller in case of cold wind solution. The 
mass out flow rate for this solution is Anq = 8.1 Mq yr~^, half 
the value which obtains for the hot wind solution. Therefore the 
cosmic ray pressure can by itself drive an outflow, with significant 
mass loss, even if the gas efficiently radiates away its thermal en- 
ergy. Note that were we to lower the adopted value of the parame- 
ter Ce, keeping all other parameters the same, one would also get 
a valid solution, with lower values for Pc, p (x Ce and B oc \/C7l. 
The velocity profile will however not change. Lowering Ce would 
also imply a lower value of q, with q oc Ce. 



5.3 Cold wind solution for small galaxies 

We now turn to the small mass galaxies. In Fig.|6]we show the cold 
wind solution (by thick lines) for a galaxy of mass 10^ A^©, adopt- 
ing the parameters as in model E except for Ce. We adjust Ce so 
that the cosmic ray pressure for the cold and hot wind solution of 
Fig. [3] are similar at the critical point. This hot wind solution is 
also shown in Fig.|6]as thin lines, for comparing with the cold wind 
solution. From the figure, it is clear that the qualitative features of 
the cold wind solution are similar for a small mass galaxy and a 
massive galaxy. The flow again starts with a small velocity near the 
centre of the galaxy potential. The final velocity is larger than that 
of hot wind solution but still related to the circular velocity with 
F/V2 = 1.25. The mass outflow rate is -iirq = 0.38 M© yr'^ 
and is again a factor 2 lower for the cold wind compared to the 
hot outflow. However the exact reduction in the mass outflow rate 
again depends on the parameters of the solution adopted as dis- 
cussed above for the case of winds from massive galaxies. 

Recall that the mass loading factor scaled as r]w oc V^^ for 
the hot wind solutions without radiative cooling. We have seen that 
the mass outflow rate for the cold wind solutions scales with the 
mass outflow rate for the hot wind. Thus we expect that the scaling 
relation between the mass loading factor rj^ and the circular veloc- 
ity, rjm oc V^^, will continue to hold for the cold wind solutions, 
albeit with a smaller absolute value of q, for any particular galaxy 
mass. In order to make a more quantitative statement, one needs to 
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solve for wind solutions with cooling included in the energy equa- 
tion. This is beyond the scope of this exploratory work. 



6 DISCUSSION AND CONCLUSIONS 

We have studied here cosmic ray driven wind solutions relevant to 
galaxies whose gravitational potential is of the NFW form. Cosmic 
rays are expected to be copiously generated in high redshift star 
forming galaxies, (and also in present day starbursts), where the 
SFR and hence the SNe rate, is likely to be much higher than those 
in our Galaxy. The CR driven winds extensively studied in the con- 
text of our Galaxy, could then also be relevant for such galaxies. 
Moreover, transonic wind solutions, with the fluid velocity being 
accelerated from subsonic to supersonic speeds exist for a 7 = 5/3 
thermal gas, only if one can have energy and/or momentimi input 
till the critical point. A purely thermally driven wind solution as- 
suming such energy input upto the critical point is given by Cheva- 
lier and Clegg (1985). On the other hand, cosmic rays naturally 
provide such energy and/or momentum input, even if the energy 
sources (i.e. SNe) are well within the critical point, or the thermal 
energy is lost due to radiative cooling. These considerations formed 
the key motivation for our work. 

We have first given a suite of such galactic free wind solutions, 
driven by both thermal and cosmic ray pressures, in an NFW grav- 
itational potential, where the physical parameters of the wind are 
fairly reasonable. The solutions start at radii smaller than the criti- 
cal point, with ro/l ^ 0.4 — 0.8, depending upon the dimension- 
less parameters. For galaxy masses in the range 10* — 10^^ M©, 
the gas and cosmic ray pressures at the critical point range re- 
spectively from about 2 — 0.5 eV cm~^, while B decreases from 
13.2 — 6.2 fiG. These values are comparable to the cosmic ray en- 
ergy density of our Galaxy and magnetic field strengths in nearby 
spirals. The density of the wind material at r = /, which ranges 
from ~ 2.1 — 0.003 cm~^, for the same range in galaxy masses, 
is also comparable to the densities expected in galaxies. Corre- 
spondingly the temperature of the wind increases from about T ~ 
1.5 X 10*K for a 10* Mq galaxy to 2.6 x 10® K, for a 10^^ M® 
galaxy. The increased density and decreasing temperature associ- 
ated with these winds for lower mass galaxies imply that gas cool- 
ing would become more important for such galaxies. 

It appears that /jG strength magnetic fields are required for 
CRs to drive a transonic wind in the NFW potential that we have 
adopted. We do find that it is possible to find cold wind solutions 
with very small magnetic fields in a point mass potential. This issue 
is of interest as Zweibel (2003) has shown that even fields several 
orders of magnitude smaller than galactic fields can accelerate and 
confine CRs. The magnetic field that is required, need not however 
be coherent over large scales, as we discussed earlier. Radio obser- 
vation of several edge on galaxies do show significant radial fields 
emanating from the galaxy, and associated with outflows (Heesen 
et al, 2009a, b; Krause, 2008). 

All these solutions importantly have an asymptotic wind ve- 
locity closely related to the circular velocity of the galaxy. Such 
a correlation seems to be indeed observed in outflows of ultralu- 
minous infrared galaxies (Martin, 2005 Fig. 7). Martin (2005) also 
finds that wind velocity is proportional to the SFR of the galaxy. 
In our models, larger mass galaxies have a larger SFR, and also a 
larger circular velocity and hence wind velocity. Thus one naturally 
expects to find such a correlation in the CR driven wind models. 

Our wind solutions also predict an anti-correlation between 
the mass loading factor of the wind and the circular velocity, with 



rjiu/fccVc ^ oc v^^l- We find that larger mass galaxies with M ~ 
IQii-iQi^ Mq, typically have 7?^ ~ 0.5-0.2, adopting / = 0.1. 
Such mass loading factors are consistent with 0.1 < rjw < 0.7 
determined by Rupke, Veilleux & Sanders (2002) for ultraluminous 
infrared galaxies at redshifts 0.04 < z < 0.27. For smaller mass 
galaxies with total mass M ~ 10^'' — 10^ Mq, we obtain a larger 
rjw ~ — 7.6, while for galaxies with M ~ 10* Mq this increases 
to r]w ~ 28. 

Note that in our wind models the wind temperature varies be- 
tween 0.6 — 2.2 X 10 ' K, for dwarf galaxies with total mass in the 
range 10^ — lO'^^M®. The wind can have larger temperatures for 
smaller values of the dimensionless parameter 'a' (by a factor upto 
4) and smaller values of the parameter T (the location of the critical 
point). Thus galactic free winds in our models can have tempera- 
tures required to produce X-ray emission only in massive galaxies 
but not in dwarf galaxies. However, several local dwarf galaxies do 
show extended X-ray emission indicative of hot gas, with tempera- 
ture of order KeV (cf. Martin, 1999; Ott et al., 2005), even though 
starbursting dwarf galaxies are only a small fraction (about 6%) of 
local low-mass population (Lee et. al. 2009). If this X-ray emission 
is related to the wind material then this hot component of the wind 
cannot be explained by the models discussed here, where the source 
region is confined to be within the sonic point. 

If one does wish to construct very hot free wind solutions 
in such dwarf galaxies, one will have to use a model Uke that of 
Chevalier and Clegg (1985), where the sources do extend at least 
upto the radius where the wind velocity becomes sonic. Moreover, 
for such solutions, the radiative losses cannot be more than about 
30% (cf. Silich, Tenorio-Tagle and Rodriguez-Gonzalez, 2004). It 
may be that both these conditions are obtained for some dwarf 
galaxies. But they need not hold in general; i.e. the sources need 
not in general extend to the sonic point and the physical conditions 
could be such that radiative cooling is important. In such cases, our 
work shows that Cosmic rays could still help to drive a wind. 

Indeed, in our cold wind solutions, we consider models where, 
due to radiative cooling, the thermal pressure can be neglected com- 
pared to the cosmic ray pressure. We find that even under these cir- 
cumstances, cosmic rays can nevertheless drive a healthy transonic 
wind. The flow starts at a small radius with a small velocity which 
become supersonic beyond a critical point, reaching an asymptotic 
value of order the circular velocity of the galaxy. We show that in 
such solutions, and depending on the model parameters, the mass 
outflow rates are smaller, by at least a factor 2, compared to the hot 
CR driven wind solutions. Nevertheless, even in this case, we still 
expect that smaller mass galaxies have a higher r^^. A large rj^ for 
small mass galaxies imply that cosmic ray driven outflows could 
provide a strong negative feedback to the star formation in dwarf 
galaxies. 

We have assumed spherical symmetry in our wind solutions. 
It will be useful to find such solutions in disk like geometry as well. 
We have also not included the effect of CR diffusion. Such diffu- 
sion could lead to a reduction of the mass outflow rate as shown by 
Dorfl (2004). Further, we have assumed the wind to be free, and not 
confined by any external pressure, say due to gas in the halo of the 
galaxy. In case such confining agents are important, one could first 
have a stage of filled hot bubble evolution, as in models described 
in our earlier work on thermally driven outflows into the halo/IGM 
(Samui, Subramanian & Srianand, 2008). And the free wind so- 
lution would become relevant in latter stages, when the confining 
medium has been evacuated out. 

The results on small mass galaxies obtained here, that their 
winds have asymptotic speeds close to their circular velocity, and 
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have a large mass loading factor, would have important implica- 
tions for models of metal enrichment of the IGM. Several groups 
have shown that outflows from such small mass galaxies are pos- 
sibly the most important contributors to the volume filling of the 
IGM with metals (cf. Madau, Ferrara and Rees, 2001; Samui, Sub- 
ramanian & Srianand, 2008). Thus it is important to revisit such 
global models in the light of the results obtained in this work. 
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APPENDIX A: THE WIND EQUATION 

Using pvr^ = g in Eq.js] we obtain, 
' P 



pv 



dv 5 d 
dr 2 dr 



GM{r) , 

= -pv 5 [v + Va) —r— 



Multiplying Eq.|2]by v, and subtracting from Eq. |Al| we get, 

dPc 



(Al) 



dP 
dr 



dr 



5p d_ /P 
Tdr \j> 



{v + va) 



dr 



(A2) 



which can be simplified to give 

dP _ 5Pdp 
dr 3 /9 dr 

The Alfven velocity is given by, 



VA = 



2 VA dPc 

3 u dr 



(A3) 



B 



(A4) 



where we assume that the magnetic field satisfies the condition 



which leads to 



Br 



1 dVA 

Va dr 



constant, 



1 dp 
Ypdr' 



Also Eq.[T]can be differentiated to get 



1 dp 1 dii 2 
p dr V dr r 



(A5) 



(A6) 



(A7) 
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Hence from Eq.|4] Eq. |A7| and Eq. |A6| we get the cosmic ray pres- 
sure gradient as, 



dr 3 (u + va) L 2 J p dr 



(A8) 



The gas pressure gradient can be obtained using Eq. |A3| and Eq. |A8 

"5P 8Pc (u + ua/2) 



dP 
dr 



dr 



3 p 9 p (?; + va) 

Substituting Eq. |A8| and Eq. |A9| in the momentum equation (Eq.[2]l, 
we get, 



Av 2 dp pGM (r) 



dr 



' dr 



where the effective sound speed c, is determined by, 

c2 = 5P _^ 4P, {v + va/2){2,v-2va) 
3 p 3 p ?,v{v + va) 



(AlO) 



(All) 



Eliminating dp/dr using Eq. |A7[ we get the wind equation given 
in the main text. 



Av 
dr 



2d 



GM{r) 



2rci 



(A12) 



APPENDIX B: HOT WIND WITH UNIFORM SOURCE 
DISTRIBUTION 

We consider here the hot wind solution in presence of energy 
sources upto certain radius R. This will allow us to estimate what 
fraction of total SNe energy ends up as the kinetic energy of the 
wind material. We assume that upto r — R, SNe explosions 
provide energy to the thermal gas and cosmic rays as well. Let 
Cth and Ccr are the luminosity per unit volume injected to the 
thermal gas and cosmic rays respectively. We take both Cth and 
Ccr as constant for r < R and zero out side. These two source 
terms have to be added to the right hand side the energy equation 
(Eq. [3]l and the cosmic ray evolution equation (Eq. |4j. Note that 
(4/3)7rP^(£th + Ccr) = (4/3)7rP^£o = Lq, where Lq is the 
total luminosity that is being pumped by the SNe to the wind mate- 
rial. Apart from energy, SNe explosions would also add mass and 
momentum to the system. Hence the modified fluid and cosmic rays 
equations can be written as (compared to Eqs.fTpl, 



J__d_ 

^2 (Jj, 



(pvr^) 



Q = constant if r < P 
if r > P 



(Bl) 



where Q is the mass input rate per unit volume to the system. Fur- 
ther, Eq.|5]for the energy exchange between cosmic rays and ther- 
mal gas remains the same. 

One can integrate Eq.|Bl|to get 



pvr 




if r < P 
if r > P. 



(B5) 



Note that, q is the same constant as in Eq.[T] Demanding the conti- 
nuity of p and V across r = P, we obtain. 



QP^ 



(B6) 



We can also add Eqs. |B3| and |B4| and integrate the resulting 
equation to get an algebraic conservation law. For r > P the source 
free equations are still valid and we will get the conservation law, 
Eq.[6] However, for r < P, we get, 



pvr 



1 2^ 5P 
-v H 

2 2 p 



+ 



QGM 
35(c) 



fvir 
C 



2 1 

4 



cr 

fvir 



In 1 + 



CT 

fvir 



+APcr^ (v + va)^^ 



(B7) 



Again demanding continuity across r = R, from Eqs. [6] and 
|B7|(and also using the relation of Eq.[B6}, we obtain. 



47rCe 



SnqGM 
d{c)R 



2 1 



~cR 



In 1 + 



cP 



(B8) 



The second term in the right hand side comes from the gravitational 
term only. Note that AttCc is the total energy flux that is passing 
through any sphere and it is also the final wind kinetic energy. So, 
the final kinetic luminosity of the wind is the total luminosity that is 
being supplied by the SNe minus the energy per unit time required 
to climb out of the gravitational potential of the galaxy. Hence the 
fraction / of Eq.[TT] entirely depends on the gravitational potential 
of the galaxy if radiative cooling is negligible. If the total luminos- 
ity is less than the required potential energy, there would not be 
any outflowing solution. This threshold luminosity would depend 
on the shape of the gravitational potential and the distribution of 
the energy sources with respect to that potential. 

In order to determine / for NEW galaxy potential and the 
uniform source distribution, we assume that Lo is higher than the 
threshold value and hence inCc is positive. We use Eqs.|l l|and|12| 
to rewrite Eq.|B8|and find 
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if r < P 
ifr>P ' 
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(B3) 



(B4) 
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S r 1)2. 
1 + - - —HO 



(B9) 



where 



HO = 2C (1 - C') In 1 + ^ + C (2C - 1) 



(BIO) 



with 



c = 



Note that Vdr ~ Vc for a nearly flat rotational curve. 

The function is shown in Fig. |B1| It is clear that J-" is of 



Cosmic ray driven outflows 15 




1 2 3 4 5 



Figure Bl. The function of Eq. [BTo] is shown as a function of It reaches 
to an asymptotic value of 1.3 as f oo. 



order unity for a wide range of ^ values. For example, if we take 
R = ry^r/12 and c = 15, we get C = 0.8 and then T = 0.94. 
Further, 5(c) — 1.84 for c = 15. We will see f rom numerical 
solution that F/ \f2 ~ 1. Putting these values in Eq. B9 we obtain 
/ = 0.08. Hence most of the total luminosity has been used to 
climb out of the gravitational potential and only about 10% ends 
up as wind kinetic energy. 



